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We propose a new shell model method, combining the 
Lanczos digonalization and extrapolation method. This 
method can give accurate shell model energy from a series 
of shell model calculations with various truncation spaces, in 
a well-controlled manner. Its feasibility is demonstrated by 
taking the fp shell calculations. 



PACS numbers: 21.60.-n, 21.60.Cs 



Quantum many-body phenomena frequently reveal 
novel and emerging nature beyond a simple extrapo- 
lation from the single-particle picture, mean-field pre- 
dictions, and exact diagonalization of Hamiltonian ma- 
trices by such as the Lanczos method. Solving such 
many-body problems has been a common challenging 
issue in elementary particle physics, nuclear structure 
physics, condensed matter physics and others. In these 
decades, diverse numerical methods, such as quantum 
Monte Carlo, stochastic diagonalization/variation, den- 
sity matrix renormalization group (DMRG), were pro- 
posed. Among them, as a method combining the Path- 
integral formalism and diagonalization/renormalization 
group, in nuclear structure physics, the Quantum Monte 
Carlo diagonalization (QMCD) method [QJ|] was pro- 
posed and was developed. In condensed matter physics, 
the Path-integral renormalization group (PIRG) method 
H was proposed. Though both share common proce- 
dure, in the latter method, an extrapolation method is 
utilized. This extrapolation method is significant by it- 
self and can be applied, independently of other parts of 
the PIRG method. In the present paper, we combine this 
extrapolation method with the Lanczos diagonalization 
and propose a new shell model method. 

The shell model is one of the fundamental frameworks 
for nuclear many-body problem. In the conventional shell 
model calculation, Hamiltonian matrix elements of the 
effective shell model interaction are calculated by a com- 
plete basis set in a given shell model space, and then, di- 
agonalization is carried out by the Lanczos method. Such 
Lanczos shell model calculations have a long history over 
a half century and have extended its realm of the ap- 
plication due to the progress of computer technologies. 
Nowadays, part of full fp shell problem can be solved, 
for instance, the yrast states of 52 Fe, the M scheme di- 
mension of which in the complete fp shell model space is 
109,954,620, was studied Q, while there still exist many 
nuclei not yet explored in the full fp shell due to their 
huge shell model dimensions. 

For such nuclei, truncation scheme is utilized for reduc- 



ing the shell model space. In the fp shell calculations, 
due to relatively large gap of spherical single particle en- 
ergy between the fj/2 orbit and others, we restrict the 
shell model space by specifying the number of excited nu- 
cleons across this shell gap, that is, ©s<t(/V/2)' 4-40 ~ s ( r ) s 
where r means the set of the fs/2, P3/2 and px/2 orbits 
and t is the maximum number of nucleons allowed to 
be excited. As an example, we consider 48 Cr within fp 
shell model space, of which shell model calculation needs 
about 2 million dimension. This was one of the state-of- 
the-art large-scale shell model calculations in mid 1990's 
II . In Fig. 1 (a) and 1(b), the ground state energies 
of the truncated space are plotted as functions of the t 
value and the corresponding M scheme dimensions, re- 
spectively. The definite behavior of the convergence is 
not seen in Fig. 1, because energy eigenvalue as a func- 
tion of the i or M scheme dimension does not have a 
well-defined scaling property. Therefore, in general, we 
cannot predict the ground state energy in the entire shell 
model space by extrapolating shell model results within 
smaller truncated spaces unless the truncation calcula- 
tions give well-converged results. Hence, it is crucial to 
find a useful extrapolation scheme to extract the correct 
results from truncation calculations. 
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FIG. 1. Ground state energies of 48 Cr as functions of the 
t values (a) and the M scheme dimensions (b). The labels t 
are added for clarifying the truncation spaces in (b). 

In order to extrapolate the shell model results within 
smaller truncated spaces into the complete shell model 
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result, a well-defined scaling property for energy eigen- 
values is needed. We define the difference SE between 
the energy eigenvalue (H) in a given truncated space and 
true energy eigenvalue {H) g , that is, SE = (H) — (H) g . 
The energy variance AE in a given truncated spaces also 



defined as, AE — 



The difference SE vanishes 



linearly as a function of the energy variance AE. This re- 
lation was utilized in Refs. j| and |(| and the summary of 
its proof can be seen in Ref . Q . An approximate ground 
state \ip) can be decomposed by the true eigenstate \ip g ) 
and rest component \ip e ) as \tp) = c \ip g ) + d \ip e ) where 
c 2 + d 2 = 1, and \ip g ) and \ip e ) are orthnormalized states. 
Up to 0{d 3 ), SE (x AE is satisfied @. Therefore, this re- 
lation shows a good scaling if is a good approximation 
of the true eigenstate. 
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FIG. 2. Extrapolation of the energy to the zero energy vari- 
ance for the 0^ and 2+ states of 48 Cr. The t values represent 
the truncation space. 

As an example, we reconsider the shell model results 
of 48 Cr by the energy variance AE as the a;-axis. In Fig. 
2, energy eigenvalues for various truncation spaces are 
plotted as a function of the energy variance for the 0^ 
and 2\ states of 48 Cr with the KB3 interaction ||. Lin- 
ear relation between SE and AE is seen from 4=1. By 
the results corresponding to the t = 1,2,3, t = 3,4,5, 
t = 4, 5, 6, we can extrapolate the ground state energy 
as -33.148 ± 0.15, -33.047 ± 0.02, -32.975 ± 0.02 MeV, 
respectively, where these error-bars are given by the \ 2 
fitting. As the exact ground state energy is -32.954 MeV, 
the extrapolation from t = 1,2,3 results is already a 
good value. Moreover, as the truncation space becomes 
larger, the error due to the extrapolation procedure be- 
comes smaller. Thus the extrapolated energy is given 
with error-bar, and its error can be well-controlled. In the 
SMMC calculation a different extrapolation method 
is utilized for overcoming the minus sign problem [ jl0| . 
The extrapolated energy of this nucleus is -32.3 ± 0.4 
MeV O]. This error comes from two reasons. One is a 
statistical uncertainties due to the Monte Carlo sampling, 
the other is an systematic error due to the extrapola- 
tion and finite temperature corrections. Compared with 
the SMMC calculation, the present method has only a 
small systematic error. Here, we comment that, extrap- 



olated energies may become lower than the exact energy 
because this process does not necessarily guarantee the 
variational principle. 
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FIG. 3. Extrapolation of the quadrupole moment (a), oc- 
cupation number of the / 7/ / 2 orbit (b) for the 2+ states and 
B(E2:2+ -> 0+) (c) to the zero energy variance. The effective 
charges are taken as e p = 1.5e and e„ = 0.5e. 

Next we consider the expectation value of various oper- 
ators by the present extrapolation method. In Fig. 3, we 
show the quadrupole moment, occupation number and 
B(E2,2+ -> 0^) for 48 Cr. As the first two values are ex- 
pectation values by the same wavefunction, the AE can 
be uniquely defined. We can extrapolate them into zero 
energy variance for quadrupole moments. In Fig. 3(a), 
unlike energy, t = 1,2 results seem to be rather poor. 
From the t = 1, 2, 3, t = 3, 4, 5 results, the extrapolated 
quadrupole moments are -24.6 ± 4.1, -29.0 ± 0.3 efm 2 , 
respectively. The latter extrapolated value reaches al- 
most exact one (-29.1 efm 2 ). On the other hand, for the 
B(E2), the initial and final states are different. Therefore, 
initial and final states have different AE's, which makes 
extrapolation procedure difficult. In this case, however, 
the shell model spaces with the same t value but differ- 
ent total magnetic quantum number have similar AE. 
Therefore, average AE works quite well in Fig. 3(c). 
The extrapolated B(E2) from t = 2, 3, 4 results is 220 ± 
7 e 2 /m 4 , which very well reproduces the exact one ( 220 
e 2 /m 4 ). 

By this example, we conclude that energy vs. energy 
variance plot is useful to get an extrapolated energy from 
a series of shell model results with various truncated 
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spaces. In addition to the energy estimates, the same 
extrapolation is also shown to be possible for quadrupole 
moments, B(E2) and other quantities. 

Moreover we comment on the present extrapolation 
from the structure of 48 Cr. Low- lying yrast states of 
48 Cr is deformed §, and for its description, the particle- 
hole excitation across the shell gap between the m orbit 
and others is essential. Therefore it is interesting that the 
extrapolation from the truncated shell model results with 
few particle-hole excitations works quite well. 



> 
S 



a 



■ 54 Fe 

-170 

°i 

■□ 

-175 


t=l o 

(a) 

t=5~^^ 
t=4 


8 
6 

□ 

: * 


1 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 L 
- (b) / / ^ =1 " 

jit t=2 

t=3 

<"t=5 



0.0005 0.001 0.0015 

Energy variance 

FIG. 4. (a) Extrapolation of the energy to the zero energy 
variance for 54 Fe with other results of the SMMC |uj (rilled 
circle with error bar), statistical truncation |l^](open circle), 
t — 7 truncated shell model (open square) and the QMCD ]14] 
(filled square), (b) Extrapolation of the total Gamow Teller 
strength to the zero energy variance for 54 Fe with other re- 
sults of the extrapolation concerning to the t value |l6[ (open 
circle), t = 7 truncated shell model (open square) and the 
SMMC ^ (filled circle with error bar). 

We proceed further large-scale shell model problem. 
The M scheme dimension of 54 Fe is about 0.35 billion 
and several recent methods |h]|1^] solved this shell model 
problem. For astrophysical interest, the Gamow Teller 
transition of this nucleus is known to be important. The 
Gamow Teller transition also offers good benchmark test 
for the present extrapolation. For the following calcula- 
tion of 54 Fe, the FPD6 interaction jl3| is taken as the 
effective interaction. Fig. 4(a) shows the ground state 
energies of t =1,2,3,4 and 5 shell model spaces. The ex- 
trapolated ground state energy from t = 3, 4, 5 is -174.49 
± 0.02 MeV. Here we compare this value to the ener- 
gies of other methods. The ground state energy of the 
t = 7 shell model spaces, of which M scheme dimension is 
91,848,462, is -174.217 MeV. The extrapolated energy of 
the SMMC calculation go| is -172.9 ± 0.6 MeV. In Ref. 
fl2fl , from the statistical point of view, the JT scheme 
bases are selected and then the Hamiltonian matrix is 
diagonalized. This energy is about -172.9 MeV O]. The 



ground state ene rgy of the QMCD method with ten J- 
projected bases [Q is -173.9 MeV. The present extrap- 
olated energy is consistent with the QMCD and t = 7 
truncated shell model energies. 

In Fig. 4(b), the total GT+ strength is plotted. The 
extrapolated total GT+ from t = 3, 4, 5 result is 4.98 ± 
0.06. The total GT+ strength of t = 7 calculation is 
5.03. In Fig.2 of Ref. [jj), the total GT+ strength is 
plotted as a function of the t value and it suggests that 
its total strength is about 5. On the other hand, in Ref. 
p0| , the extrapolated value of the SMMC calculation is 
4.32 ± 0.24. Our extrapolated value is consistent with 
the former two values. 

The 56 Ni is a key nucleus for understanding the shell 
structure around the fp shell region because it has a dou- 
bly magic structure, while its closed shell is suggested to 
be rather soft (IJitJ]- Moreover, recently two deformed 
bands are discovered ]Is[| . This nucleus is quite interest- 
ing but the M scheme dimension is quite large ( about 
1.1 x 10 9 ). Then this nucleus was a target of a state- 
of-the-art shell model calculation in the end of 1990's 
|p|, p^|p8| . Therefore it is also a good touchstone for the 
present extrapolation method. 
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FIG. 5. Extrapolation of the energies of the 0J~, 2+ , 2 J and 
2g" states to the zero energy variance for 56 Ni. For each state, 
the results of the t = 2 ~ 6 truncated shell model spaces are 

slotted. In the left part, the corresponding QMCD results 

151 are shown. 



In Fig. 5, the energies of the Of, 2j~, 2J 



and 2g states 



with t — 2 ~ 6 shell model spaces are plotted. In princi- 
ple, the present extrapolation method can be applied to 
the excited states with the same quantum numbers with- 
out any difficulty. Here we test its numerical feasibility. 
The extrapolated ground state energy from t = 4 ~ 6 is 
-203.26 ± 0.06 MeV. In Ref. Q, the ground state energy 
of the QMCD calculation is -203.152 MeV. The ground 
state energy of the t = 8 truncated space is -203.178 MeV 
p0f, of which M scheme dimension is 255,478,309. The 
extrapolated ground state energy is consistent with the 
results of the QMCD calculation and the t = 8 truncated 
shell model calculation. The extrapolated 2f energy is 
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-200.46 ± 0.05 MeV from t = 4 ~ 6 results, which is 
also consistent with the QMCD result and the t = 8 
truncated shell model result (-200.075 MeV) @. 

Next we investigate the excited 2 + energies. In the left 
part of the Fig. 5, the QMCD energies |L5| are shown. 
These 2 + states are characterized by the quadrupole mo- 
ments and the occupation numbers in the orbit. The 
quadrupole moments are 21.6, -41.9 and 9.5 e/m 2 for the 
2j, 2 2 and 2% states, respectively. The occupation num- 
bers are 12.9, 11.6 and 12.1. In Ref. the 2\ state is 
concluded to be a prolate deformed state. As the trunca- 
tion scheme is based on the spherical single particle en- 
ergies, the deformed 2^~ state is considered to be difficult 
to describe within small t truncated spaces. For t < 6, 
the states with negative quadrupole moment appear as 
the third 2 + state. By plotting the shell model ener- 
gies for the t = 2 ~ 6 truncated spaces, we can find the 
linear relations for 2y , 2^ and 2^~ energies, respectively, 
as shown in Fig. 5. Therefore, for excited states, the 
present extrapolation method can be applied. Moreover 
Fig. 5 shows crossing between the 2^ and 2^ energies 
as a function of the energy variance, which means that, 
the present extrapolation method can correctly handle 
these excited states with different nature. The extrap- 
olated energy of the 2^~ state from t = 3 ~ 6 results is 
-198.35 ± 0.35 MeV and extrapolated quadrupole mo- 
ment and occupation number are -28.2 ± 8.5 efm 2 and 
11.5 ± 0.6, respectively, which are consistent with the 
QMCD values. For more presice description of this pro- 
late deformed state, t = 7 results are needed, while it is 
quite surprising that such a deformed excited state can 
be described by the present extrapolation method based 
on the small truncated shell model spaces. 

Finally we discuss a computational aspect. For the 
present extrapolation method, in addition to the usual 
Lanczos procedure, the evaluation of the matrix element 
(H 2 ) is required. To obtain this matrix element, we oper- 
ate the Hamiltonian H to the l^) of the t truncated space 
in the [t + 2) truncated space, then we can evaluate the 
norm of the H \ip). Therefore, the implementation of this 
extrapolation method is quite easy for a usual shell model 
program. 

In summary, we have proposed a new shell model 
method based on the Lanczos diagonalization and an ex- 
trapolation method. By taking the state-of-the-art fp 
shell calculations in the last decade, we have shown its 
numerical feasibility. The present extrapolation method 
can give accurate energies, quadrupole moments and 
other quantities from a series of shell model results with 
smaller truncated shell model spaces. 

The present extrapolation method was successfully ap- 
plied to condensed matter problems as a part of the 
PIRG method |22|. Therefore, it is intriguing to examine 
whether this extrapolation method is useful also in the 
QMCD method for the nuclear shell model. For instance, 
we expect more precise energy estimates from the QMCD 
calculations with smaller number of bases. Moreover, it 



might become quite useful in other methods of nuclear 
structure physics. Such researches are in progress. 

We thank Professor N. Yoshinaga for reading the 
manuscript and useful comments. This work was 
supported from 'Research for the Future Progam' by 
Japan Society for the Promotion of Science under grant 
RFTF97P01103, and was supported in part by Grant-in- 
Aid for Scientific Research (A)(2) (10304019), and one for 
Specially Promoted Research (13002001) from the Min- 
istry of Education, Science and Culture. 



[1] M. Honma, T. Mizusaki and T. Otsuka, Phys. Rev. Lett. 

75, 1284 (1995); Phys. Rev. Lett. 77, 3315 (1996). 
[2] T. Otsuka, M. Honma and T. Mizusaki, Phys. Rev. Lett. 

81, 1588 (1998). 
[3] M. Imada and T. Kashima, J. Phys. Soc. Jpn. 69 

2723 (2000). 

[4] C. A. Ur, et al., Phys. Rev. C58, 3163 (1998). 
[5] E. Caurier et al, P hys. Rev. Lett. 75 2466 (1995). 
[6] S. Sorella: preprint, |cond-mat/0009149 
[7] T. Kashima, M. Imada, J. Phys. Soc. Jpn. 70 2287 
(2001). 

[8] A. PovesandA. P. Zuker, Phys. Repts. 70, 235 (1981). 
[9] S. E. Koonin, D. J. Dean, and K. Langanke, Phys. Repts. 
577, 1 (1996). 

[10] Y. Alhassid, D.J. Dean, S.E. Koonin, G. Lang, W.E. 

Ormand, Phys. Rev. Lett. 72 613 (1994). 
[11] K. Langanke, D.J. Dean, P.B. Radha, Y. Alhassid, 

S.E. Koonin, Phys. Rev. C52, 718 (1995). 
[12] M. Horoi, B. A. Brown, V. Zelevinsky, Phys. Rev. C50 

R2274 (1994). 

[13] W. A. Richter, M. G. van der Merwe, R. E. Julies and 

B. A. Brown, Nucl. Phys. A523 325 (1991). 
[14] T. Mizusaki, T. Otsuka, M. Honma and B. A. Brown, 

Phys. Rev. C63 044306 (2001). 
[15] T. Mizusaki, T. Otsuka, Y. Utsuno, M. Honma and 

T. Sebe, Phys. Rev. C59, 1846R (1999). 
[16] E. Caurier, G. Martinez-Pinedo, A. Poves, and A. P. 

Zuker, Phys. Rev. C52, R1736 (1995). 
[17] G. Kraus, et at, Phys. Rev. Lett. 73, 1773 (1994). 
[18] D. Rudolph, et at, Phys. Rev. Lett. 82, 3763 (1999). 
[19] T. Otsuka, M. Honma, T. Mizusaki, N. Shimizu, and 

Y. Utsuno, Prog. Part. Nucl. Phys., 47 319 (2001). 
[20] These are the largest shell model calculations in our shell 

model code MSHELL §l|. 
[21] T. Mizusaki, RIKEN Accel. Prog. Rep. 33 14 (2000). 
[22] For example see, T. Kashima and M. Imada, J. Phys. 

Soc. Jpn., 71, 3052 (2001); H. Morita, S. Watanabe and 

M. Imada, unpublished. 



4 



